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O ■ ABSTRACT 

The metallicity in galaxy clusters is expected to originate from the stars 
in galaxies, with a population dominated by high mass stars likely being the 
<^ \ most important stellar component, especially in rich clusters. We examine the 

fr^ I relationship between the metallicity and the prominence of galaxies as measured 



in 



by the star to baryon ratio, M^:/Mf,ary Counter to expectations, we rule out 
a metallicity that is proportional to M*/Mbary, where the best fit has the gas 

O I phase metallicity decreasing with M^/Mhary, or the metallicity of the gas plus 

the stars being independent of M^/Mhary This implies that the population of 
stars responsible for the metals is largely proportional to the total baryonic mass 

'k>( I of the cluster, not to the galaxy mass within the cluster. If generally applicable, 

H ' most of the heavy elements in the universe were not produced within galaxies. 

Subject headings: galaxies; clusters: general — galaxies: clusters: intracluster 
medium — X-rays: galaxies: clusters 



1. Introduction 

While nucleosynthesis is undoubtedly responsible for the production of the elements, 
there is less certainty regarding the stellar populations involved, as well as the time and lo- 
cation of elemental production. Historically, these issues have been studied within individual 
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galaxies, an ideal approach if they are nearly close d box systems. H owever, studies find that 
galaxies have lost the majority of their baryons (JMcGaughl l2005l ) and some undetermined 
amount of metals. For example, a galaxy like the Milk y Way is missing 75% of its o riginal 
baryon content while for M33, the deficit is about 90% ([Anderson fc Bregman 1120101 ). 



Galaxy clusters provide another useful system for studyi ng metallicity evolut i on, a s 



moderate and rich clusters retain most of their baryons (e.g., iBohringer fc Wernerll2009l ). 
Baryons lost from gal ax:ies are captured in the hot cluster medium, which is effectively 
studied in X-rays (e.g., IWerner et al.l 120081 ). The first order of business is to assess whether 



the stars in galaxies have produced the observed metallicity. Simple estim ates indicate a 



problem in that there are too few stars to accou nt for the heavy element mass (jArimoto et al. 



19971 : iPortinari et al.ll2004t iLoewensteinl 120061 ) . We estimate that the metallicity of stars in 
galaxies, averaged over a standard luminosity function, is about 0.5 solar. Galaxies comprise 
about 5-10% of the baryons in rich clusters, so if the metals in these stars were spread 
through the cluster, we expec t the hot gas raetalli city to be about 5% solar, but the typical 
metallicity is 40% solar (e.g., ISnowden et al.ll2008l ). about an order of magnitude higher. 



This issue was examined by several authors ( iPortinari et al.l |2004J : ILoewensteinl |2006| ) , 
who did careful calculations of metallicity production from the var ious types of super - 
novae, using acosmic star formation history (review of techniques by lBorgani et al.l 120081 ). 
LoewensteinI ( 120061 ) concludes that with any of the standard initial mass functions, the metals 
produced by the stars in galaxies cannot produce the cluster metals. He posits the existence 
of a population of mainly high-mass stars to pro duce the cluster metallicity and he inves- 
tigates the time e volution. IPortinari et al.l (120041 ) reach a si milar conclusion and adopt the 



IMF proposed by lArimoto et al.l (119971 ). which is top-heavy ( Fabian et al.l l2008l use a top 



heavy IMF and a normal IMF). A dissenting view is given by ISivanandam et al.l ( 120091 ) who 
find that for some clusters, the metals might have originated in galaxi es . However , s evera l 
of t heir adopted values ar e relatively extreme compared to the work of ILoewensteinl (120061 ) 
and IPortinari et al.l (|2004[ ) . 



In a related paper I Anderson fc Bregman I (|2010| ). we argued that field spirals are missing 
their baryons due to the presence of an early stellar population dominated by high-mass 
stars, which heats and enriches the gas before most of it collapsed into proto-galaxies. In 
this picture, which is driven by simple observational constraints, galaxies are baryon-poor 
because the gas never fell into galactic dark matter potential wells. Preheating (and pre- 
enrichment) by this early population would need to occur near the onset of galaxy formation, 
which occurs before cluster formation. Therefore, this picture should apply to systems that 
eventually become galaxy clusters. 



It is natural to expect that the mass of the high-mass weighted early population of stars 
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is proportional to the mass of a galaxy, so the collective mass of this population in a cluster 
should be proportional to the total stellar mass seen in galaxies today. Translating this into a 
metallicity prediction, the metallicity of the cluster should be proportional to the ratio of the 
stellar mass to the total baryon mass, M^,/M^,ary■ This is the same expectation if the metals 
were produced and expelled from the cluster galaxies. We show that this basic expectation 
is excluded by the data, potentially leading to a profound change in our understanding of 
metallicity production in clusters and in the universe. 



2. Object Selection and Data Products 

For this study, we need to measure four parameters of galaxy clusters: the gaseous mass; 
the stellar mass; the total gravitating mass; and the metallicity. From these quantities, we 
can assess the total baryon fraction of clusters, along with the ratio of the stellar mass 
to the total baryonic mass, and compare it to the metallicity. These quantities are most 
easily measured for low redshift clusters, which are the objects used. There are a variety of 
measurements of the properties of low-redshift clusters, but only a limited number measure 
the mass to large radii (rsoo) and are close enough to measure the stellar mass content from 



sources such as the 2 MASS data set (JKochanek et al.ll2003l : iLin et al.ll2004J : IVikhlinin et al. 



2009 : iDai et al.ll2010l ). We discussed many of these objects in a separate work ( iDai et al. 



2010l ). from which we select clusters that have retained most of their baryons, so have 
similar total baryon fractions (Table 1); these clusters appear to be relaxed. The gaseous 
masses are determined by fitting a functional form to the surface brightness distribution. The 
functional form begins with a /3-model, which can be modified so that it steepens at large radii 
plus another component in the center (IVikhlinin et al.ll2009l ). Temperatures are determined 
from the data using spectral fits (i.e., APEC models), and this temperature information 
is important i n determining the t otal g ravitating mass of the system. Metallicities are 
normalized to lAnders &: Grevessd (119891 ) where the relative abundances are held at these 
solar ratios. In determining a metallicity, the most prominent lines are from Fe, which is the 
most important element statistically. Henceforth, we use [Fe/H] as a proxy for metallicity, 
where the solar elemental ratios are assumed. A ACDM concordance cosmology is assumed 
with Ho = 73 km s"^ Mpc-\ fi^ = 0.26, and ^a = 0.74. 

Metallicity has been obtained as a function of radius for many of our systems by us 
and others, from which one can i nfer an emission- weighted metallicity or a mass-weighted 
metallicity (lAnderson et al.l l2009l and references therein). The difference is generally not 
large, so we use the emission-weighted metallicity, as this is the only quantity available for 
clusters that are more gas-poor. The primary complication in making this measurement is 
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accounting for the presence of cooling cores, which can bias the emission-weighted metaUicity 
upwards by introdu cing additional flux from th e enriched core, usually coincident with a 



central galaxy (e.g.. lRasmussen &: Ponmarul2009l ). We therefore present two metallicities for 



each cluster, one emission-weighted over all of the X-ray emission, and one that excludes 
emission from the core of the cluster, as described below. 

Correctly accounting for a cool core requires careful deprojection of the emission, but 
we lack the data for such an analysis in all cases, and so per form a simpler correction. We 
examined the universa l temperature and abundance profiles of iRasmussen &: Ponmaru (120071 ) 
for galaxy groups and iBaldi et al.l (120071 ) for hot (T > 6 keV) galaxy clusters. Both groups 
found temperature profiles that slowly rise towards the core, then decline more sharply 
starting at i? ~ 45T^/^ kpc, where T is the emission-weighted temperature of the ICM 
in keV. Most of the objects in our sample have temperatures within the range of values 
considered in these papers, so the same profiles should also apply. We therefore attempt to 
exclude emission from within a projected radius of 45T^/^ kpc in each cluster to produce 
an emission- weighted iron abundance without the bias of a bri ght core (T is in keV units). 
This is straightforward for the clusters with measurements from lSnowden et al.l (120081 ). since 
they provide measurements of the temperature, metaUicity, and flux in a series of annuli 
for each cluster, which allows us to verify that a cooling core exists and then to remove the 
inner annuli and recompute the flux-weighted metaUicity. We perform a similar analysis on 
3C442A, for which we have measurements of temperature and metaUicity as a function of 
ra dius (Sun 2009, p rivate communication), and for RGH80, using the single-temperature fits 
of IXue et al.l ( 120041 ). For these two systems, we estimate the expected flux in each annulus 
by assuming the ICM gas density follows a /^-model with /3 = 0.65. In Abell 1275 the 



meta llicities are from our analysis of the Chandra data. For A160, and A2462 ( jJetha et al. 



20051 ). we have only global measurements of the metaUicity, so we estimate the effect of a 
cooling core. Based on the other clusters in our sample, the the exclusion of a bright core 
reduces the iron abundance by 25%, so we use this correction. 

The final quantity is the stellar mass, which is obtained by identifying the galaxies and 
measuring their magnitudes to a radius of r2oo; which includes nearly all of the stellar light. 
The data set for the galaxies used here is the 2MASS data base, so there is a magnitude 
limit to galaxy identification. This leads to sampling only part of the galaxy luminosity 
function, so a correction is applied to account for galaxy i ncompleteness and Poisson bias 



for a simple halo occu pation model, as previously discussed ( JKochanek et al.ll2003l : iLin et al. 



2004J : iDai et al.l 120071). The uncerta inties assigned are conservative in that they are larger 
than those given in lLin et al.l ( 120041 ). For some of the analysis, we include the metals in the 
stars, so we assumed a solar metaUicity for that correction as applied to the total stellar 
mass. 
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This procedure does not include a cont ribution from intracluster light (ICL). Measures 
of the intracluster light were determined by lKrick fc BernsteinI (120071 ) for 10 galaxy clusters, 
where their fractional ICL results are consistent with a single value of 11.1 ± 1.3% (x^ min- 
imization). For individual clusters, the B and r weighted ICL values range from 4.8 ± 3.1% 
for Abell 2556 to 21.3 ± 6.6% for Abell 4059 and there is no obvious correlation of the ICL 
percentage w ith cluster propertie s. This measurement agrees with the value from theoretical 
modeling by iRudick et al.l ( l2006l ) an d with the frac t ion of host-less sup ernovae in c l usters , 
which is quoted as 20% (+20/-12) bv lCal-Yam et al.l J2003h and ~20% bv lSand et all J2008h . 



Gonzalez et al. 



(120071 ) ■ but they combine the ICL 



Another study that includes ICL is that of 
and the light from the brightest central galaxy, so this does not provide an independent 
measure of the ICL. They imply that the ICL may increase toward poorer clusters. In the 
following analysis, we discuss the effects of ICL that comprise 10-20% of the stars, as well 
as an increase of the ICL toward poorer clusters. 



3. Metallicity as a Function of Baryon Fraction 

With the set of clusters for which we have good data, we examined the relationship 
between the metallicity of the gas and the fraction of the baryons in the form of stars (Fig 
1). A proportional relationship between these quantities is the nominal expectation, but the 
best fit, Z = 2.44 yi^/lsAi^ary^ has a badly unacceptable y^ = 729 for 10 dof, so it is ruled 
out. Upon removing RGH 80, the object with the highest value of M^,/Mf,ary, the fit is still 
ruled out, as x^ = 413 for 9 dof. The data are characterized by a metallicity that is slowly 
declining as a function of cluster mass, while the importance of the stellar mass increases 
toward lower cluster mass. 

We consider another model to explain the data, one where part of the metals are associ- 
ated with stars B(M^./Mbarj/), and the other where some fraction of the metals has an origin 
independent of the current visible stars, Zq. This second component can be understood as an 
early population component (e.g., dominated by high mass stars) that is not proportional to 
the galaxy mass. This is expressed as a linear equation, Z = Zq + B (M^^/Mbary)- We made 
fits to two sets of metallicities, a single emission-weighted metallicity for the cluster and a 
metallicity with the core excluded. The linear fits in the two cases are similar, but the fit is 
better for the metallicities without the core, so we concentrate on those results (Fig. lb, 2). 
For the case without the metal contribution from the stars, the best-fit is a negative slope, 
Z = 0.37 - 0.70{MjMbary), with errors on Zq and B of 0.02 and 0.23; x^ = 18-5 for 10 dof, 
which is acceptable. For the case with the stellar metal contibution included, the best fit 
has a fiat slope, Z = 0.37 + 0.10{MjMbary), with errors on Zq and B of 0.02 and 0.22; x^ = 
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14.4 for 10 dof. The second fit is consistent with a constant metalhcity, but the fit without 
the stellar metalliity is inconsistent with a constant metallicity solution. 

The inclusion of ICL either does not change our result or make it more restrictive. When 
calculating t he stellar compon ent, we used a K band M/L of 0.95. Had we used the IMF 
suggested by IChabrierl (120031 ) . the M/L would be about 0.72, resulting in a stellar mass 
that is 24% lower. On the oth e r han d, inclusion of ICL would increase the stellar mass by 
11%, using iKrick &: Bernstein! ( 120071 ). Either of these effect lead to minor scaling changes 
in yi^/yibaryi but do not change the results. If the ICL increases significantly toward poorer 
clusters, which preferentially occur in the high M^,/Mf,arj/ part of the diagram, the maximum 
allowable slope for the gas metallicity is lowered as is the fraction of metals than can be 
associated with galaxies. 

This general result, that the metallicity is constant or decreases as the stellar fraction 
increases, has been known for over a decade. An increase in the ratio of the stellar mass 



to the gas mass with decreasing cluster temperature was first demonstrated by iDavid et al. 



( 1l990l ). A r nodest decrease of clus ter metallicity with decreasing cluster temperature was 
discussed by iFukazawa et al.l (119981 ). Our work confirm these trends but uses more uniform 
stellar metallicities, gas masses, and total gravitating masses, all measured to large radii. For 
our conclusions to be invalid, a large body of work by others would need to be dramatically 
wrong. 

Although we used systems in which the baryons are mostly retained by the clusters, the 
situation is not dif ferent if extended to galaxy gr oups in general. In their recent survey of 15 
groups of galaxies, iRasmussen fc PonmanI ( l2009l ) find that the metallicity of the gas beyond 
the central galaxy is lower than in clusters, a trend that is consistent with our result (Fig. 
lb). 

The implications of this result are remarkable in that it suggests that not only are most of 
the metals in the richest clusters produced by the an IMF dominated by high mass stars, but 
that the amount of mass in this component is not proportional to the stellar mass in galaxies. 
The need to produce most of t he metals in an IMF that is heavily weighted on the high-mas s 
end was previously discussed (JArimoto et al.lll997l : iLoewensteinl l2006l : IPortinari et al.l 120041 ). 
This inference is supported by the elemental ratios within rich c lusters, which favor meta l 
production by an ensemble of high-mass metal-poor progenitors ( jBaumgartner et al.ll2005l ). 
The mass of this stellar population would appear to be roughly proportional to the mass 
of gas in a cluster. This conclusion would explain the mean cluster metallicity, but not the 
metallicity gradient within a cluster. For that to be explained with this stellar component, 
the frequency of such SNe would need to be biased toward higher density regions. This is 
similar to the arguments that were put forward for biased galaxy formation, so this might 
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be considered biased star formation with a high-mass IMF. 

Clusters of galaxies were used in this study because they represent a "closed box" piece of 
the universe. The situation in clusters is likely to be representative of the universe as a whole, 
in which case, most of the heavy elements in the universe were produced outside of galaxies 
by an early population of supernovae. If correct, this represents a fundamental change in 
our understanding for the origin of the elements. This moderately enriched material (0.1-0.4 
solar) would fall into galaxies, producing a disk with few low-metallicity stars (a previously 
proposed solution to the G-dwarf problem). The increase of the metallicity within galaxies 
is of course due to star formation, mass loss and supernovae within the galaxy. 



In addition to the above issues, some additional evidence supports our model. lEhlert Sz Ulmer 

0.89 and from the lack of 



(120091 ) analyzed cluster metallicity profiles from z = 0.14 to z 
evolution, they conclud e that the metallicity distribution was established at high redshift. 
David fc NulsenI ( l2008l ) note that Fe is more extended than the light of stars, which could 
imply that the galaxies and the metals ar e not closely coupled (they sugge st an alternative 
model where AGN heating is responsible). ISommer-Larsen fc Fynbd (120081 ) find that at z ~ 
3, less than 20-25% of the oxygen is associated with galaxies, which is also consistent with 
most of the metals being produced outside of galaxies. 

If there was a predominantly high-mass population of stars outside galaxies, the initial 
mass function of these stars could be constrained by the ICL, provided one could separate 
the ICL from galaxy-galaxy interactions with that from the high-mass stars. The present 
values of the ICL must provide a limit to the slope of that IMF, but future work may 
be able to constrain it more directly. Another constraint on the na ture of the metal pro- 
duction has been the va r iation of the elemental ratio s with radius (JDupke fc White! l2000l : 
Baumgartner et al.ll2005t iRasmussen &: PonmanI 120091 ). In groups, the Si/Fe ratio increases 
with distanc e away from the central gala xy, indicating that Type la are less important at 
larger radii (IRasmussen &: PonmanI |2009| ) . There is significant room for improvement in 
these types of studies, by using more elements, and especially for the hot clusters. 

In obtaining this result, only a dozen clusters were used, a shortcoming that can be 
rectified in future work. Rich clusters are well-represented and the number of such objects 
can be increased with straightforward archival work of X-ray data; these tend to have low 
values of M^^/Mbary The greater weakness is the number of galaxy clusters with high values of 
M^/Mhary Thcsc tcud to be of lower mass, relatively gas-poor and cooler systems, so the X- 
ray emission is weaker, making metallicity determinations from X-ray data more challenging. 
Additional X-ray data will be needed to make progress here. Finally, more accurate values 
of M* can be determined with deep imaging, as typically only two dozen galaxies define L* 
in a cluster. While these are the the brightest galaxies, the fainter end of the luminosity 



function should be determined in each case. 
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Fig. 1. — Metallicity in the gas as a function of the stellar fraction of baryons, M^,/Mf,ary, 
for a sample of 11 clusters with similar total baryon fractions, measured relative to the 
cosmological value. For the top panel, the metallicity is measured for the whole cluster while 
for the lower panel, the metallicity is measured with the cooling core excluded. The steep 
dotted line, which is ruled out, is the model predicted if the metallicity is proportional to 
the galaxies and passes through the rich clusters (the four clusters on the left). The red 
line is the best-fit linear relationship and the two curved green lines are the 95% confidence 
bounds. The results imply that very little of the metals in rich clusters (low M^,/Mhary) is 
due to galaxies. 
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Fig. 2. — This is the same as the above, except the metalhcity of the stars is added as well 
(solar metallicity assumed for the stars). The metallicity contribution from the gas is with 
the cooling core excluded. This data set also rules out a model where the metallicity is 
proportional to the stellar fraction of baryons, M*/Mfcarj/- 



Table 1. Galaxy Cluster Mass and Metallicity Contents 



Cluster 


Alt Name 


Redshift 


T 
keV 


M. 
M© 


Mtot(r5oo) 
Mq 


b/rac 


cr 


M./M^ary 


(J 


T^tot 


CF 


^nc 


(J 


Abell 133 




0.0566 


4.14 


6.03E+12 


3.13E+14 


0.112 


0.007 


0.072 


0.017 


0.57 


0.09 


0.45 


0.09 


Abell 478 




0.0881 


7.94 


1.20E+13 


7.57E+14 


0.156 


0.014 


0.042 


0.012 


0.40 


0.02 


0.36 


0.02 


Abell 1795 




0.0625 


6.12 


4.57E+12 


5.95E+14 


0.132 


0.007 


0.024 


0.006 


0.45 


0.02 


0.37 


0.02 


Abell 1991 




0.0587 


2.61 


3.72E+12 


1.21E+14 


0.141 


0.010 


0.091 


0.024 


0.52 


0.06 


0.37 


0.06 


3C442A 




0.0263 


1.61 


1.19E+12 


3.90E+13 


0.100 


0.006 


0.129 


0.025 


0.30 


0.07 


0.28 


0.07 


NGC 5098 


RGH80 


0.0362 


1.05 


2.20E+12 


2.00E+13 


0.190 


0.021 


0.251 


0.064 


0.20 


0.04 


0.15 


0.04 


Abell 160 




0.0447 


1.99 


3.45E+12 


7.40E+13 


0.132 


0.010 


0.150 


0.032 


0.40 


0.10 


0.30 


0.10 


Abell SllOl 


Scrsic 159/03 


0.058 


2.69 


1.51E+12 


1.41E+14 


0.147 


0.025 


0.030 


0.010 


0.40 


0.02 


0.30 


0.02 


Abell 1275 




0.0603 


1.63 


3.72E+12 


6.90E+13 


0.144 


0.030 


0.159 


0.055 


0.40 


0.04 


0.34 


0.05 


Abell 2462 




0.0733 


2.62 


2.04E+12 


8.80E+13 


0.135 


0.012 


0.072 


0.026 


0.24 


0.05 


0.18 


0.05 


Abell 2029 




0.0773 


8.47 


1.15E+13 


8.64E+14 


0.157 


0.011 


0.047 


0.003 


0.50 


0.06 


0.38 


0.05 



Note. — Only M-toti^Mo) is listed at rsoo. AH other quantities are corrected to r200. M, is the stellar mass, Mjot is the total gravitating mass, 
hfj-ac is the baryon fraction, M*/Mi,u^y is the stellar to total baryon ratio, Ztot is the emission-weighted metallicity for the whole cluster, and Z„c is 
the metallicity without the central cooling region. 
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